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Abstract 

Ginzburg-Landau expansion for superconductors with anisotropic s— and 
d— wave pairing is derived in the presence of anisotropic normal impurity 
scattering which makes d— pairing state more stable under disordering. It is 
demonstrated that the slope of the upper critical field \dH C 2/dT\T c in super- 
conductors with d— wave pairing has nonlinear dependence on disorder, i.e. for 
the low anisotropic scattering rate it drops rather fast with concentration of 
normal impurities, but as anisotropy of scattering increases it features initial 
nonlinear growth, approaches a maximum and drops again, vanishing at the 
critical impurity concentration. In superconductors with anisotropic s— wave 
pairing \dH C 2/ dT\x c grows, approaching the known asymptotic behavior, char- 
acteristic of usual isotropic case irrespective to the scattering anisotropy. 

PACS numbers: 74.20.Fg, 74.20.De 
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I. INTRODUCTION 



The main problem of the present day physics of high-temperature superconductors is 
the determination of the nature and type of Cooper pairing. A majority of experiments 
and theoretical models |J suggest the realization in these systems of anisotropic pairing 
with (i :r 2_ y 2-symmetry with zeroes of the gap function at the Fermi surface. At the same 
time there exist some theoretical and experimental evidences m supporting the so-called 
anisotropic s— wave pairing. In this latest case there again appear zeroes (with no change of 
sign) or rather deep minima of the gap function in the same directions in the Brillouin zone 
as in the case of d— wave pairing. It was shown that controlled disordering (introduction 
of normal impurities) can be an effective method of experimental discrimination between 
different types of anisotropic pairing. In our previous paper it was shown that one can use 
measurements of the slope of the upper critical field \dH C 2/ (IT\t c for the same purpose, i.e. 
for the case of d— wave pairing the slope drops rather fast with disorder, while for the case 
of anisotropic s— wave pairing it grows with disorder as for the isotropic case. 

Recently an interesting theoretical model was proposed |7|] taking into account the ef- 
fects of anisotropic (momentum-dependent) impurity scattering. It was shown that for large 
enough anisotropic "d— wave" scattering the usual pair-breaking effect of normal impurities 
(described by Abrikosov-Gorkov dependence for isotropic superconductor with magnetic im- 
purities) is rather strongly suppressed. This allows, in principle at least, to overcome one 
of the main problems in the physics of high-temperature superconductors — the contra- 
diction between the d— wave picture of pairing in these systems and their relative stability 
towards disordering M . There exist also some other explanations of such a stability (cf . an 
explanation proposed in our paper ||), however the simplicity of the suggested model @ 
is rather attractive and calls for further theoretical study of superconductors with exotic 
pairing with the account of anisotropic impurity scattering. The present paper is a natural 
continuation of our previous study, including the effects of anisotropic impurity scattering. 
It is demonstrated that anisotropic impurity scattering leads to significant anomalies in the 
dependence of the slope of the upper critical field on disorder (impurity concentration). As 
in the previous paper || our analysis is based on the microscopic derivation of Ginzburg- 
Landau expansion in impure system. 



II. GINZBURG-LANDAU EXPANSION 

Following Refs. 0|| , we analyze two-dimensional electronic system with isotropic Fermi 
surface and separable pairing potential of the form: 

V(p,p') = V(<j>,<t>') = -Ve(<f>)e(<!>'), (1) 

where is a polar angle, determining the electronic momentum direction in the conducting 
plane, and e(4>) is given by the following model dependence 

e(6) = I ^ C ° S ^ ( d_ WaVe )' (2) 
1 \/2\cos (2(f)) | ( anisotropic s— wave ). 

The pairing constant V is as usual different from zero in some region of the width of 2lj c 
around the Fermi level (u; c -is some characteristic frequency of the quanta, responsible for 
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the pairing interaction). In this case the superconducting gap (order parameter) takes the 
form: 

A(p) = A(0) = Ae(0), (3) 

and positions of its zeroes for s and d cases just coincide. 

We examine a superconductor containing "normal" (nonmagnetic) impurities, which are 
chaotically distributed in space with concentration p. Following Ref. |F] we consider the 
square of the scattering amplitude of the impurity in the following form: 

IWp,p')I 2 = IW^OI 2 = |v o r + M7(0)/(^), ( 4 ) 

where Vq is isotropic point-like scattering amplitude, V\ is anisotropic scattering ampli- 
tude, and /(</>) is angular-dependent model function (0 is a polar angle mentioned above) 
which describes the type of anisotropic scattering. We consider the scattering to be essen- 
tially isotropic and impose the following constraints 0: 

N 2 <N 2 ; </>=o ; </ 2 >=i, (5) 

where < ... > denotes the average value over the momentum direction on the Fermi surface 
(i.e. over the 0-angle). Accordingly, the second part in Eq.(f4|) represents the deviations from 
the isotropic scattering. 

The normal and anomalous temperature Green's functions in such a superconductor are 

G(u, p) = , (6) 

& + % + \A(p)\*' 

F( W ,P) A * (P) 



^+^+iA( P )r 

where u = (2n + l)vrT, £— is electronic energy with respect to the Fermi level, 

u{p)=uj + ip J ^^\V imp (p - p')\ 2 G{uj,p'), (7) 
A(p) = A(p) +p J ^- 2 \V imp (p-p')\ 2 F(uj,p'). 
To determine the transition temperature we can limit ourselves to Eqs. (|7|) linearized over 

u = co + ip ^- 1 dt / o 27r #{|v | 2 + N70/W)} (8) 

a = a + p^M f di r d<p{w + \vi\ 2 mm) ^72- 

Z7i J Jo L J uj z + 1, 2 

The critical temperature T c is determined by the linearized gap-equation: 

d P' „, * A (P' 



A: 



a(p)--^E/^(p.p')^ o» 



Following standard procedure from Eqs.(|8]),(|9|) we obtain the following general equation for 
the critical temperature T c : 



In 



T 

- L c 

Tea 



< e > 2 + < ef > 2 -1 



<ef >' 



* - - 





7o \ 




2nT c J 


7o 


fl-- 


2uT c 


V 7o 



(10) 



where T c0 - is the transition temperature in the absence of impurities, *f?(x)- is the usual 
digamma function, 7 = 7rpVg 2 iV(0) and 71 = npV^N(0) - correspondingly the isotropic and 
anisotropic impurity scattering rates, < ef > 2 describes the "overlap" between the functions 
e(p) and /(p). 

For the sake of simplicity we take the function /(p) in the form analogous to that in 
Eq.(§: 



/(p) = f(<f>) = V2cos(2ct>), 



This corresponds to the maximum overlap for d-case. More general treatment one could find 
in Ref. f?J. Now the renormalized Eqs.(H) can be written as follows: 



to = uj + i 



7T 

• 7o 



UJ 



A = A + i— I d£— 



uj 2 + £ 2 
A 



+ i-^cos(2<p) [d£ ! d<p'cos(2<p')- 

7l 2 J J UJ 



■ 7l 



7T Z 



+ i^cos{2<f)) / d£ / d(p' cos {2<f)')- 



u 2 + e 2 ' 

A 



(12) 



From here we obtain the well-known expression for the renormalized frequency in both cases: 

uj = uj + j signuj. (13) 

In the case of d— wave pairing the gap symmetry in the presence of impurities is not 
changed: 



A = A 





Co 




\uj 




7i 



(14) 



In the case of s— wave pairing the gap is shifted by a constant, which does not depend on 
and 71: 



A = A + A 



2%/2 7o 



7T M 



Finally T c -equation for superconductor with d— wave pairing is written as: 



T, 



c0 



Jl\ 7o \ 
7o/2tttJ- 



For superconductor with anisotropic s— wave pairing T c -equation is written as: 



' n 777- 
T, 



c/0 



l ~V 2 



2nT r 



ml 



(15) 



(16) 



(17) 



Note that anisotropic scattering rate dependence drops from Eq.fll7|). 

The appropriate dependencies of T^o/T^o) are shown in Fig.l, for the case of d— wave 
pairing with different values of the normalized anisotropic scattering rate 71/70- In the case 
of s— wave pairing the transition temperature T c is slightly suppressed with the growth of 
7o/T c o. In the case of d— wave pairing T c suppression is rather fast for small values of 71, 
but the critical value of 7oc/^co leading to the complete destruction of superconducting state 
rapidly increases with the growth of the anisotropic scattering rate 71/70- 

The gap function as usual can be used as an order parameter in the Ginzburg-Landau 
expansion for the free-energy density. We consider its amplitude A(T) to be a slowly varying 
function of the spatial coordinates. Accordingly in momentum space we get the Fourier- 
component of the order parameter: 

A(0,g) = A g (T)e(0). (18) 

The Ginzburg-Landau expansion for the free energy density difference between super- 
conducting and normal state in the region of small q up to terms quadratic over A takes the 
form: 

F s -F n = A\A g \ 2 + q 2 C\A g \ 2 (19) 

and is determined by the diagrams of loop-expansion for the free-energy of electrons moving 
in the field of superconducting order parameter fluctuations with some small vector q, shown 
in Fig. 2. Diagrams (c) and (d) are to be subtracted, so that the coefficient A becomes zero 
at the transition point T = T c . 

Some details on calculating r pp / and Ginzburg-Landau coefficients for the case of d— wave 
pairing can be found in the Appendices A and B. Note, that for the d— wave superconductors 
the contribution of diagrams Fig.2(b,d) actually vanishes up to terms of the order of q 2 , 
if we don't take into account an anisotropy of impurity scattering. In the case of s— wave 
superconductor all calculations are analogous, we only note that in such dependence 
on anisotropic scattering rate is absent. 

Finally we can express GL-coefficients in the following form: 

A = A K A ; C = C K C , (20) 
where Aq and Cq are just the usual GL-coefficients for the case of isotropic s-wave pairing 
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T - T c „ ,^ 7 C(3)^ 



A = N(0)-^; C = N(0)^^, (21) 

where vf, iV(0)-are electron velocity and normal density of states at the Fermi level, and 
everything specific to our model is contained in dimensionless combinations Ka and Kc- In 
the absence of impurities for both models we obtain: K\ = 1, Kq = 3/2. For the impure 
system we get: 

(A) d-w&ve pairing: 

Ka = r d 4 r du A + (22) 
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(23) 



(24) 



A, 



28tt 2 C(3) V2 2ttT c 
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24vr 2 



T 2 



7C(3)7o 2 (^ 2 -8)7 2 VT c0 ; 7C(3) 7 o 



(25) 



The appropriate dependencies of dimensionless coefficients on disorder parameter 7o/T c0 in 
the case of d— wave pairing and for different values of the normalized anisotropic scattering 
rate 71/70 are shown in Figs. 3,4. 



III. UPPER CRITICAL FIELD 



GL-coefficients A and C, as usual, define the temperature dependence of the upper 



critical magnetic field close to T c : \TL 



H, 



c2 



2^ 2 (T) 



~2itC 



(26) 



where 0o = cn/e - is magnetic flux quantum, £(T) - is temperature dependent coherence 
length. Now we can easily find the "slope" of the temperature dependence of H C 2(T) near 
T c , i.e. the temperature derivative: 



dH, 



c2 



dT 



247r0o m K A 



7((3)v 2 F C K C 



(27) 



In the case of the usual s-wave superconductivity anisotropic scattering does not influence 
the behavior of the slope of the upper critical field. The appropriate dependencies of dimen- 
sionless parameter h = \dH c2 /dT\T c /\dH C 2/dT\T c o on disorder 7o/T c0 in the case of d— wave 
pairing - for different values of the normalized anisotropic scattering rate 71/70 are shown 
in Fig. 5. In the case of anisotropic s-wave pairing the slope as usual grows with disorder 
and in the limit of strong disorder 70 3> T c it crosses over to the usual linear dependence 



described by the well-known Gorkov's expression [12]: 



N(0) 



dH, 



c2 



dT 



8e 2 J 

7T 



(28) 



where a = N(0)e 2 v 2 7 /3 , jo - is electron conductivity in the normal state, which is charac- 
teristic of the impure superconductors with isotropic s-wave pairing. It means that strong 
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disordering suppresses gap anisotropy and we obtain a usual limit of the impure supercon- 
ductor. 

For the case of <i-wave pairing the slope of H c2 drops to zero on the scale 70 ~ T c0 for 
the small values of the rate 71/70- For the values of anisotropic scatterong rate on the 
interval 0.5 < 71/70 < 0.6 the behavior of the slope is qualitatively changed: h smoothly 
but nonlinearly increases with the growth of 7o/T c0 , crosses over the maximum and then 
has a sharp drop. The interval where the slope grows extends as 71 approaches 70. In 
our opinion these sharp anomalies in dependence of the slope of the upper critical field 
on disorder can be used for determining the pairing type and possible role of anisotropic 
impurity scattering in "unusual" superconductors. Unfortunately, in case of high-T c oxides 
the situation is complicated by the known nonlinearity of temperature dependence of H C 2, 
which is observed in rather wide region close to T c and and also by some ambiguity in 
experimental methods to measure H c2 . 

This work was partly supported by the grant No. 96-02-16065 of the Russian Foundation 
of Fundamental Research. It was performed under the project IX. 1 of the State Program 
"Statistical physics" as well as the project No. 96-051 of the State Program on HTSC of the 
Russian Ministry of Science. 
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APPENDIX A: VERTEX PART r PP , IN "LADDER" APPROXIMATION. 

The Bethe-Salpeter equation for the vertex part takes the form: 

r pp > = U(p, p') + E ^(P, p")G R (p")G A (p")T p „ p ,, (Al) 
p" 

where U(p, p')-is irreducible vertex function. We take U(p,p') in the following form ("lad- 
der" approximation): 

^(p,pO = p^ 2 + pv; 2 /(p)/(p'). (A2) 

Then Eq.(Al) can be written as: 

r pP ' = P V 2 + P V?f{p)f{p') + pV^{p') + plf/(p)$(p') (A3) 

where 

^(p') = E^(p")^(p")r P ^, (A4) 
P " 

$(p , ) = E/(p ,, )^( P ")G A (p")r P ^- 
p" 

From Eq.(A3) one can obtain a self-consistent set of equations for ^(p') and $(p'): 

/ *(p') = P v 2 h + pv; 2 /(p')/ 2 + pV 2 i^(p>) + pv; 2 / 2 $(p'), 

\ $( P ') = pv'h + pif/( P ')/ 3 + pK, 2 / 2 ^( p ') + pv; 2 /3$(p'), [Ab) 

where 

A=E^(P)G A ( P ), 
p 

/ 2 = E/(P)^(P)G A ( P ), (A6) 
p 

/ 3 = E/ 2 (p)^(p)^(p). 
p 

Solving system (A5), one can find the appropriate expressions for ^(p') and $(p') and hence 
the expression for the vertex part: 

r p^ 2 (l - pV?I 3 + pVffjpVi) + pV 1 2 (f(p)f(p')(l - pVjh) + P v 2 f(j>)h) 

(i - pK, 2 / 1 )(i - P v?i 3 ) - pvipvfit [A{) 
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APPENDIX B: GINZBURG-LANDAU COEFFICIENTS. 



We can easily see that the contribution of the diagram Fig. 2 (a) is 

T 



(2k) 



AlY, I d P 2cos 2 (2<f ) )G^p + )G^ 

/ • T J 



,^ 2 N(0)7ru|T cv , 1 

+ e ~ q 
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a; 



The contribution of the diagram Fig. 2(b) is 
T 



(27T) 



A£/ rfp2co S 2 (20)G.(p)G_(p) = - A 2 T c iV(0) £ / 



£ 2 + £ 2 ' 

The contribution of the diagram with diffusion propagator Fig. 2(b) is 

-T^^v / 2co S (20)G^(p + )G^p_)r pp ,v / 2co S (20')G R (p' + )G^(p'_). 



pp 



Taking into account (A6) and (A7) we get from here 
-Tiy(0)7r 7l ]T 



1 v 2 F (2\u\ 


-7ik 2 




-7i) 8 


CO 


3(|a 


l"7i) 2 



(Bl) 



(B2) 



(B3) 



(B4) 



Note that in the absence of anisotropic scattering for the case of d— pairing the contribution 
of diagrams (c) actually vanishes up to terms of the order of q 2 . 

In the same way we get the appropriate contribution of the diagram (d) 



-7W(0)7T 7l £- 



w\(\u>\ -7i)" 



(B5) 



Finally we get the expression for F s — F n and so Ginzburg-Landau coefficients cited in the 
main body of the paper. 
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FIGURE CAPTIONS 



Fig.l. Critical temperature T c as a function of the normalized isotropic scattering rate 
7o/T c0 . The dashed curve represents the s-wave pairing case, the solid curves represent the 
(i-wave pairing case for different values of the normalized anisotropic scattering rate 71/70: 

l—7i/7o =0.0; 2—0.3; 3—0.5; 4—0.6; 5—0.7; 6—0.8; 7—0.9; 8—0.95. 

Fig. 2. Diagrammatic representation of Ginzburg-Landau expansion. Electronic lines are 
"dressed" by impurity scattering. T is the impurity vertex calculated in "ladder" approxi- 
mation. Diagrams (c) and (d) are calculated with q = and T = T c . 

Fig. 3. Dependence of dimensionless coefficient K A /K A0 on disorder parameter 7o/T c0 . 
The dashed curve represents the s-wave pairing case, the solid curves represent the <i-wave 
pairing case for different values of the normalized anisotropic scattering rate 71/70: 

l—7i/7o =0.0; 2—0.4; 3—0.6; 4—0.7; 5—0.8; 6—0.9; 7—0.95. 

Fig. 4. Dependence of dimensionless coefficient K c /K C o on disorder parameter 7o/T c0 . 
The dashed curve represents the s-wave pairing case, the solid curves represent the <i-wave 
pairing case for different values of the normalized anisotropic scattering rate 71/70: 

l—7i/7o =0.0; 2—0.4; 3—0.6; 4—0.7; 5—0.8; 6—0.9; 7—0.95. 

Fig. 5. Dependence of normalized slope of the upper critical field h , 

on disorder parameter 7o/T c0 . The dashed curve represents the s-wave pairing case, the solid 
curves represent the ci-wave pairing case for different values of the normalized anisotropic 
scattering rate 71/70: 

l—7i/7o =0.0; 2—0.4; 3—0.5; 4—0.6; 5—0.7; 6—0.8; 7—0.9; 8—0.95. 
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